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ABSTRACT 

Substructures, expected in cold dark matter haloes, have been proposed to explain 
the anomalous flux ratios in gravitational lenses. About 25% of lenses in the Cosmic 
Lens All-Sky Survey (CLASS) appear to have luminous satellites within ^ 5 kpc 
of the main lensing galaxies, which are usually at redshift z ~ 0.2 — 1. In this work we 
use the Millennium Simulation combined with galaxy catalogues from semi-analytical 
techniques to study the predicted frequency of such satellites in simulated haloes. The 
fraction of haloes that host bright satellites within the (projected) central regions is 
similar for red and blue hosts and is found to increase as a function of host halo 
mass and redshift. Specifically, at z = 1, about 11% of galaxy-sized haloes (with 
masses between 10^^h~^ Mq and 10^^h~^ M©) host bright satellite galaxies within 
a projected radius of 5 kpc. This fraction increases to about 17% (25%) if we 
consider bright (all) satellites of only group-sized haloes (with masses between 10^^h~^ 
Mq and lO^'^h'^ M©). These results are roughly consistent with the fraction (~ 25%) 
of CLASS lensing galaxies observed to host luminous satellites. At z = 0, only ~ 3% 
of galaxy-sized haloes host bright satellite galaxies. The fraction rises to ~ 6% (10%) 
if we consider bright (all) satellites of only group-sized haloes at z = 0. However, most 
of the satellites found in the inner regions are 'orphan' galaxies where the dark matter 
haloes have been completely stripped. Thus the agreement crucially depends on the 
true survival rate of these 'orphan' galaxies. We also discuss the effects of numerical 
resolution and cosmologies on our results. 
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1 INTRODUCTION 

In the hierarchical scenario of structure formation, struc- 
tures in the Universe are assumed to have grown from tiny 
quantum fluctuations (generated during an inflationary pe- 
riod) through gravitational instability. Due to the shape of 
the power spectrum, structures form hierarchically - larger 
structures form via accretion and merging of smaller struc- 
tures. Dense cores of the smaller structures often survive 
the merging process and manifest as subhaloes in the pri- 
mary haloes. If substantial star formation occurs in these 
subhaloes (or their progenitors), then they will appear as 
satellite galaxies. 

In the Milky Way, hundreds of subhaloes are 
predicted, startiiig fro m earlier semi-analytical studies 
(IKaufFmann et al ]l 19931). to more recent high-resolution sim- 
ulation s (iKlvpin et al.fil999l: iMoore et al.1 IT999I : iGao et al.1 



a new population of satellites has been discovered in th e 
Sloan Digital Sky Survey data (e.g. iBelokurov et al.ll2007l ). 
It should be noted though that these satellites are compact 
and, in general, much fainter than the previously known 
ones, thus it is likely that even this new population of 
satellite galaxies cannot completely remo ve the discrepancy 
between simulations and observations tMadau et al. 200^). 
It is possible that many subhaloes are dark due to ineffi- 
cient star formation, for example, due to its suppression by 



cient star rormation, tor example, due to its suppression by 
the UV-background radiation (e.g.lDoroshkevich et al.lll967l : 
ICouchman fc Reesll 19861 : lEfstathioulll992l ). 



l2004al lbl: iDiemand et al.1 l2007bh ." A few years ago, there 



were only a dozen or so satellites known, far fewer than 
the predicted number of subhaloes. However, very recently. 
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Such dark substructure can, potentially, be detected 
through several means, for example through gamma-ray 
radiation due to annihilations of dark matter particles 
(iStoehr et al] l2003l : iDiemand et al.1 [2007a ). Gravitational 
lensing is, in principle, another way to detect dark (and 
lumin ous) substructure. Flux anom alies (|Mao fc Schneideij 
Il998l ). astro metric perturbations llChen et all 20o3) and 
time-delays (iKeeton fc Moustakasl l2C)08h can be used to infer 
the presence of substructure in strong gravitational lenses. 
The results of these studies are so far inconclusive (e.g. 
iKochanek fc Dalall [20o3 : iMao et al.ll2004l ). If all substruc- 
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ture is equally efficient in affecting the flux ratios, then it is 
clear that there is more than sufficient mass in subhaloes 
to explain the flux anomalies. Unfortunately, most sub- 
haloes are in the outer part of the galaxy halo, which means 
they will have relatively little impact on the flux anoma- 
lies occurring in the central parts of lensing galaxies. Curi- 
ously, as emphasised by Schneider (2007 , private communi- 
cation ), 3 of the 6 radio lenses studied bv lKochanek fc Dalall 
()2004l ) exhibit luminous satellite galaxies close to the pri- 
mary lensing galaxy, namely MG0414-I-0534, B1608-I-656 
and B2045-f265. A question naturally arises: are such lumi- 
nous satellite galaxies expected this frequently in the current 
structure formation theory? 

The lensing cross-section is dominated by elliptical 
galaxies, thus for lensing applications it is important to di- 
vide galaxies into different types, and see whether the sub- 
halo populations are different. Furthermore, we explore in 
more detail the evolution of satellite galaxies as a function 
of redshift. If the evolution is slow, we can more conveniently 
use studies of nearby galaxies to infer the properties of the 
luminous satellite population of galaxies at intermediate red- 
shift (between 0.5 and 1, where most lensing galaxies lie). 
These are the two specific aspects of the subhalo popula- 
tion that we will address in this paper. For this purpose, we 
will use the largest cosmological simulation combined with 
semi-analytical catalogues to select haloes and study their 
satellite populations. We compare our results t o the CLASS 
survey l|Browne et al.ll2003l : iM^ers et al.ll2003l ). 

The plan of the paper is as follows. In Section [2j we 
describe the Millennium Simulation and the semi-analytical 
galaxy catalogue we use. Our main results are presented in 
Section [3l and we finish with a summary in Section [l] 



2 NUMERICAL SIMULATION DATA 
2.1 Millennium Simulation 

The Millennium Simulation, run by the Virgo ConsortiurrQ, 
follows the evolution of 2160'' particles within a comoving 
box of length 500 Mpc with a force softening length 
of 5 kpc (where the Hubble constant Hq = lOOh 

km s^^ Mpc^^). A ACDM cosmology is assumed, using 
parameters consistent with the results obtained from the 
first-y ear Wilkinson Micro wave Anisotropy Probe ( WMAP) 
data (jSpergel et all |2003| ): = 0.25, = 0.75, h = 
0.73, n = 1 and as — 0.9. Haloes are identified using the 
SUBFIND algorithm to detect self bound groups of at least 
20 particles. This implies a minimum mass for a detected 
halo of 1.7 X 10^" Mq. For a fuU description of the 
simulation the reader is referred to lSpringel et al.l l|2005H ^. 



2.2 Semi-analytic galaxy catalogues 

There are three publicly available galaxy catalogues that 
have been created using the Millenn ium Simulat i on to 
trace the underlyi n g dark matter (jBower et al.l l2006l : 
iBertone et al. I l2007l : jPe Lucia fc BlaizotI ^OOtI T We have 



used the Lucia fc BlaizotI (|2007l ) galaxy catalogue 
for this analy si s. It is base d on the mo d el des cribed in 
ISpringel et al] Hooi) and ICroton et al] ||2006|) and is 
simil ar to the techniques em ployed by Springel et alj (1200 ll ) 
and De L ucia et alj (|2004| ) . T he catalogues created by 
iDe Lucia fc BlaizotI (|2007l ) and iBertone et al.l (|2007l ) are 
based on similar semi-analytic models, differing only in their 
implementation of galactic feedback. We expect that the 
results we present here would not be significantly affected 
by the differences in the models. The third model, produced 
by iBowcr et al] (|2006l ). does not differentiate between 
'orphan' galaxies and galaxies that are still associated with 
their dark matter haloes. 

For each halo identified in the simulation, a central 
galaxy is 'created' with a mass fraction in baryons of 17% 
(corresponding to the global ratio, Qb/^m, as measured 
by WMAP first-year data). Initially the 'created' galaxy 
has no stellar mass, no cold gas and zero luminosity. The 
attributed baryons are in the form of diffuse gas with 
primordial composition. The semi-analytic model uses 
merger trees taken from the Millennium Simulation to 
describe the evolution of haloes that host the galaxies. The 
formation and evolution of the galaxies is then followed by 
implementing simple physical prescriptions for the baryonic 
physics, such as gas cooling, star formation and feedback 
processes (including AGN feedback). These processes 
depend on the properties of the host halo. When two haloes 
merge, the galaxy associated with the larger halo remains 
the central galaxy while the galaxy attributed to the lower 
mass progenitor becomes a satellite. Satellite galaxies are 
stripped of their hot gas and have no new supply of cool gas. 
They are allowed to form stars until their cool gas reservoir 
is exhausted. Subhaloes are followed, after merging with a 
larger system, until the dark matter subhalo is completely 
disrupted by tidal forces. These tidally stripped 'orphan' 
galaxies are then assumed to follow the position of the most 
bound particle in the subhalo before it was disrupted, until 
it merges with the central galaxy on the dynamical friction 
timescale. More details on t he formation and e volution of 
the galaxies can be found in ICroton" et all l|2006l ). 

We note that, as the minimum mass of a resolved 
halo is 1.7 X lO^^'/i"^ M0 and the creation of a galaxy relies 
on the existence of a halo, the number of low mass galaxies 
may well be underestimated. This is unlikely to affect our 
results as the low mass satellites will have lower circular 
velocities and are likely to be faint. This is discussed in 
more detail in Section [3.31 



1 www.virgo.dur.ac.uk 

^ http:/ /www. mpa-garching.mpg.de/galform/virgo/millcnnium 



2.3 Galaxy sample 

We study the satellite population within massive galaxy- 
sized haloes (haloes with a virial mas^ between 
10^^ ^"^Mq and lO" Mq) and group-sized haloes (with 
10"/i~^Mo < Myir < 10" Mo). We do not consider 



^ The mass enclosed within the virial radius. (We define the virial 
radius to be the radius of a sphere, centred on the most bound 
particle, within which the average density is 200 times the critical 
density). 
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cluster-sized haloes as none of the CLASS lenses are found in 
such environments. The number of galaxy-sized host haloes 
considered in this analysis is around 3 x 10^, varying little 
between z = and z = 1. For group-sized haloes, the num- 
ber is aroun d 3 x lO'^ at z = decreasing to ~ 2 x 10* at 

2 = 1. As in lSales et alj l|2007l ). we have imposed a bright- 
ness cutoff of M_H < —20.5 on our central (host) galaxies to 
ensure that they have a reasonable chance of hosting de- 
tectable satellites. In any case, faint central galaxies have 
small lensing cross-sections, and will have little effect on the 
statistics (see Section [3. ip . We consider all galaxies (within 
the virial radius of their host) with Mr < —17 as luminous 
satellites. This corresponds, approximately, to a 100 particle 
halo - the morphologic al resolution limit of the simulation 
fsee lCroton et a0l2006l ). This cut would not exclude any of 
the observed luminous satellites within the CLASS sample. 
Furthermore, to aid in direct comparison with observation, 
bright satellites are required to have R-band luminosities be- 
tween 1% and 50% of that of their host. For completeness, 
we also consider fainter satellites by dropping the _R-band 
magnitude and lower luminosity ratio cuts on our satellite 
sample. 

Since the lensing cross-section is dominated by mas- 
sive, red elliptical galaxies, we wish to divide our sample 
by galaxy type and determine whether there is a signifi- 
cant difference in the subhalo population of red and blue 
galaxies. To do this, we divide our according to the B — V 
colour of the host galaxy. The galaxy populations are bi- 
modal as a function of colour, wi th a well-defined r ed se- 
quence and a blue cloud (see e.g. ICroton et al ] |2006l : their 
Fig. 9). We adopt a. B — V colour cut of 0.8 at z = 0.0, and 
use B — V values of 0.70 and 0.65 as colour cuts ai z = 0.5 
and z = 1.0 respectively. We have found that 67%, 52%, 
37% of our galaxy-sized haloes (and 97%, 97%, 94% of our 
group-sized haloes) are associated with red central galaxies 
at z = 0.0, 0.5 and 1.0 respectively. 

For nearly all of the observed lensed systems, with which 
we wish to compare, the projected (physical) separation be- 
tween the main lensing galaxy and the luminous satellite 
galaxy is about 1 arcsecond, corresponding to 4.2 kpc 
at 2 = 0.5 and 5.5 kpc at z = 1. To explore whether 
current simulations produce enough satellite galaxies in the 
inner region of a galaxy halo to explain fiux anomalies, we 
have counted all satellite galaxies (satisfying the restrictions 
outlined above) within a 5 kpc projected region from 
the centre of the host galaxy (defined as the position of 
the most bound particle) , and will refer to these galaxies as 
bright central substructures. We also explore the effect of 
increasing this to a 10 kpc projected region. 

3 RESULTS 

3.1 Theoretical predictions 

In Fig. [l] we show the fraction of galaxy-sized 
(10^^ ft"^ M0 < Myir < 10" Mq) haloes with satellite 
galaxies, satisfying our magnitude cuts, within the central 
region (projected) of the halo. The fraction is plotted as a 
function of the host galaxy's I-band luminosity (for ease of 
comparison with the CLASS data). The fraction shown is 
the mean value, when averaged over 3 independent projec- 
tions. The uncertainty corresponds to the Poisson scatter 



within each bin. The top row shows the fraction of haloes 
with bright satellites within the central 5 kpc (pro- 
jected) while the bottom row shows the fraction of haloes 
which contain any (dark or bright) substructure within the 
same region. In these plots the red population is depicted 
using filled circles, while the blue population is shown using 
squares. The three columns show the fraction of haloes con- 
taining substructure for the three different redshifts we have 
considered (left: 0.0, middle: 0.5 and right: 1.0). Note that 
the most luminous blue hosts are not necessarily the most 
massive haloes but are likely to have undergone recent star 
formation. 

We find that at z = 0, about 3 — 4% of all of our galaxy- 
sized haloes appear to have bright satellite galaxies within 
5 kpc (projected) of the centre of the host. While this 
fraction is similar in both types of galaxies, it is found to in- 
crease with redshift (rising to around 11 — 12% at z = 1). 
Extending the central region to 10 kpc (not shown) we 
find that about 10% of our galaxy-sized haloes host bright 
central substructure at z = 0, and that this increases to 
almost 27% at z = 1. If we consider all substructure within 
the projected central region, not restricting the search to 
'observable' satellite galaxies, the fractions do not change 
significantly. This is illustrated in the bottom panels, where 
we have dropped our lower magnitude and luminosity ra- 
tio cuts. The fraction of haloes containing any substructure 
within 5 kpc (projected) of the centre of the host in- 
creases only moderately from ~ 3% to about 5% at z — 0. 
This conclusion remains valid all the way up to 2; = 1, where 
the fraction increases from ~ 11% to about 15%. 

Since some of the lensing galaxies reside in 
groups (see Section 13. 2[) . we explicitly check the frac- 
tion of haloes with bright satellites within a pro- 
jected 5 kpc region in group-sized haloes (with 
10"/i"^ M0 < M■^^r < lO^^h-^ Mq). The results are 
shown in Fig. [5] The fraction of group-sized haloes with 
bright central substructure is higher than in galaxy-sized 
haloes, increasing to ~ 6% at z = and to ~ 16% at z = 1 
for red galaxies, which dominate the lensing cross-sections; 
the fraction is slightly higher for blue galaxies. 

The lensing c ross-section is roughly proportional to a'* 
(e.g. [Turner et ai][l984 ). where cr is the velocity dispersion 
of the system and, from the Faber-Ja ckson relation, L oc a** 
for ellipticals (|Faber fc Jacksonlll976l ) , thus to compare with 
observations, the fraction should be weighted by luminosity. 
The luminosity weighted fraction of hosts with bright central 
substructure is given in Table [T] At z = 0, the fraction is of 
the order of 3% for galaxy-sized haloes, increasing to about 
6% if we consider only group-sized haloes. At 2: = 1, the 
fraction for galaxy-sized haloes is about 11%, rising to 
17% for group-sized haloes, still slightly below (but possibly 
consistent with) the observed fraction of galaxies with bright 
companions (see Section [3. 2p . 

Of the systems found to host bright central substruc- 
ture, most have only one bright central satellite. Only 2% 
(3%) of galaxy- (group-) sized hosts with bright central sub- 
structure host more than one bright central satellite at 
z = 0. At 2 = 1, 5% (7%) of galaxy- (group-) sized hosts 
have multiple bright central satellites. The largest number 
of bright central satellites found within any one system is 4. 

The force softening of the simulation is 5 kpc (in 
comoving coordinates); within this region, resolution effects 
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Table 1 . Luminosity weighted fraction of galaxy-sized hosts with 
bright substructure within the central 5 kpc (projected). 

Numbers in brackets correspond to values for group-sized haloes. 







Redshift 






0.0 


0.5 


1.0 


All hosts 


3 (6)% 


7 (11)% 


11 (17)% 


Red hosts 


3 (6)% 


6 (11)% 


11 (16)% 


Blue hosts 


4 (7)% 


7 (16)% 


11 (24)% 



Table 2. Fraction of projected central satellites within the 3D 
central region of galaxy-sized hosts. Numbers in brackets corre- 
spond to values for group-sized haloes. 



''3d < force softening 


z 


= 0.0 Bright 


32 (24) % 


z 


= 0.0 Total 


26 (17) % 


z 


= 0.5 Bright 


17 (11) % 


z 


= 0.5 Total 


14 (8)% 


z 


= 1.0 Bright 


11 (6)% 


z 


= 1.0 Total 


9 (5) % 



may be significant. For this reason, we explicitly check the 
fraction of the projected subhaloes within the 3D central 
region. The fraction depends strongly on redshift. About 
32% of the luminous satellites found within the projected 
central 5 kpc are found within the 3D central region in 
galaxy-sized haloes at 2 = (for groups this decreases to ~ 
24%). At 2: = 1 this fraction is 11%, dropping to ~ 6% in 
group-sized haloes (see Table All of the satellites within 
the 3D central region are 'orphan' galaxies (see Section |3]3]). 



3.2 Substructure in CLASS lenses 

We use the CLASS survey as the primary observa- 
tional data. This s urvey discovered 2 2 new gravitationa l 
lenses in the radio JBrowne et all l2003l : [m vers et al.ll2003h . 
iKochanek fc Dalall (|2004^ used 6 lenses from this survey to 
study the amount of substructure in lensing galaxies. In this 
work, we have used the whole survey to gather statistics; 
5 of the 22 CLASS lenses have luminous satellite galaxies 
within 5 kpc of the main lensing galaxy: B1608+656, 
B2045+265, MG 0414+0534, B1127-f385 and B13594-154. 

For B1608-I-656, the main lensing galaxy (Gl) is at red- 
shift zi — 0.63. In addition, there is a faint galaxy, G2, 
about 0.73 arcsec away, which is 1.8 magnit udes fainter than 
Gl b oth in the HST F160W and F814W jKoopmans et al.1 
l2003l ) filters. There are also four groups along the line of 
sight, including one at the redshift of the lensing galaxy Gl 
(Fassnacht ct al. 2006). If G2 is at the same redshift as Gl, 
then the projected se paration is 3.4 kp c. 

For B2045-I-265, iMcKean et"al] (|2007l ) found a galaxy, 
G2, about 0.66 arcsec away from the main lensing galaxy 
Gl (at redshift 0.867), which is between 3.6 and 4.5 mag- 
nitudes fainter than the main lensing galaxy Gl in F814W 
and F160W. The photometric redshift of G2 is consistent 
with that of Gl, although it is also consistent with being at 



redshift 4 — 5. If the redshifts of Gl and G2 are the same, 
then the projected separation is 3.5 kpc. 

For the q uadru pole lens MG 0414-1-0534, 
ISchechter fc Moord 1 19931 ) found a fainter companion 
that is about 1 arcsec away from the main lensing galaxy 
(which is at redshift 0.96). The object 'X' is about 2.6-2.44 
magnitudes fainter than the main lensing galaxy in the 
HST images of F160W and F814W. If the object 'X' is at 
redshift 0.96, then the projected separation is 5.5 kpc. 

For B1127-I-385, there are al so two lensing galaxies, Gl 
and G2 (|Koopmans et all Il999l ). The fainter one, G2, is 
about 1 magnitude fainter than Gl in both F814 (/ = 22.5 
for Gl) and F555 {V = 24 A for Gl). The separation be- 
tween these two galaxies is about 0.6 arcsec. If the lensing 
galaxy's redshift is between 0.5 and 1, then the projected 
separation is about 2.5 — 3.3 kpc. 

B1359-I-154 is a six image system produced by a small 
group of galaxies. Three primary lens galaxies lie on the 
vertices of a triangle separated by 0.7 arcsec at z ~ 1 (cor- 
responding to a projected separation of ~ 3.9 kpc), 
with magnit udes in I of 22.68 ± 0.28 and 23.69 ± 0.24 and 
23.70 ± 0.33 (|Rusin et al.ll200ll ). 

In the left panel of Fig.[3]we show the luminosities of the 
observed lenses found by the CLASS survey. The patterned 
histogram shows 16 of the CLASS lenses for which we have 
redshifts and I-band magnitudes (taken from the CASTLES 
websit^B)- Over-plotted (solid histogram) are the 4 CLASS 
lenses which have been shown to host luminous satellites; for 
B1127+385, its luminosity is unknown due to the uncertain 
lens redshift. 

The right panel of Fig.[3]shows the redshift distribution 
of the CLASS lenses. All of the lenses with luminous satel- 
lites have redshifts higher than the median value of ~ 0.6. 
About 75% of the lenses with z > 0.8 have luminous satel- 
lites. We caution that the six remaining lenses of the CLASS 
sample (with unknown redshifts) may be, on average, at 
higher z. By ignoring these lenses, the redshift distribution 
may be somewhat skewed. This could mean that the proba- 
bility of a high redshift lens hosting a luminous satellite may 
not be as high as implied. 

Fig. |4] shows the difference in magnitude between the 
host and satellite galaxy versus the projected separation of 
the satellite galaxy from the host. We show a random se- 
lection of our group-sized haloes with 'dark' substructure 
(crosses) and bright substructure (circles) from our z — 1 
sample. The 5 CLASS lenses found to have luminous satel- 
lite galaxies are plotted with solid circles; for B1127-I-385, 
the horizontal bar shows the range of separations when the 
lens redshift is varied from 0.5 to 1. Selection effects may be 
complicated and have not been taken into account in this 
study. While it will be difficult to observe satellites with 
large magnitude differences at small separations, we find 
that there are also few simulated satellites at very small 
separations. (The increase in number with separation is due 
to the larger area considered). As illustrated with the his- 
tograms in Fig. |4] we find that our sample of host galaxies 
and their luminous satellites is comparable to the observed 
galaxies in the (small) CLASS sample. 



* http://www.cfa.harvard.edu/castles/ l lKochanek et all) 
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10 10.5 11 11.5 12 10 10.5 11 11.5 12 10 10.5 II 11.5 12 



logLi(L) logLi(L) logL,(L) 

Figure 1. The percentage of galaxy-sized haloes (10^^h~^ Mq < M^ir < W^'^h~^ Mq) which have substructure within the inner region 
of the halo, as a function of luminosity. The top row shows the fraction of haloes with bright satellites within the central 5 kpc 
(projected). The bottom row shows the fraction of haloes which contain any substructure within the central 5 kpc (projected) region. 
In these plots the red population is depicted using filled circles, while the blue population is shown using squares. The three columns 
show the fraction of haloes containing substructure for the three different redshifts we have considered (left: 0.0, middle: 0.5 and right: 
1.0). The Poisson scatter is shown. Note that the most luminous blue hosts are not necessarily the most massive haloes but are likely to 
have undergone recent star formation. 
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Figure 2. As in Fig.lTlbut for group-sized haloes {10^^h~^ TsAq < M„ir < 10^*/i~^ Mq). The top row shows the fraction of haloes with 
bright satellites within the central 5 kpc (projected). The bottom row shows the fraction of haloes which contain any substructure 
within the central 5 kpc (projected) region. In these plots the red population is depicted using filled circles, while the blue population 
is shown using squares. The three columns show the fraction of haloes containing substructure for the three different redshifts wo have 
considered (left: 0.0, middle: 0.5 and right: 1.0). The Poisson scatter is shown. 
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log Li(L_) Redshift 

Figure 3. Histogram showing tlic distribution of the 16 CLASS lenses with available redshifts and I-band magnitudes, while the solid 
histogram shows the distribution of the 4 CLASS lenses with lumi nous satellites (B1127+385, is not shown due to its uncertain lens 
redshift). K-correction values have been taken from iPogeiantil l(l993), interpolated using a polynomial fit. 



Table 3. Percentage of bright satellite galaxies without a surviving dark matter subhalo within the virial radius and within the central 
5 h^^ kpc (projected) region for galaxy-sized hosts. Numbers in brackets correspond to values for group-sized haloes. 







Redshift 






0.0 


0.5 


1.0 


Satellites without DM subhalo 


74 (68)% 


81 (76)% 


87 (83)% 


Projected central satellites without DM subhalo 


98 (98)% 


99 (98)% 


99 (99)% 



3.3 Resolution effects 
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Figure 4. Difference in I-band magnitude between the host and 
satellite galaxy versus the projected separation (in kpc) of 
the satellite galaxy from the host. We show a random selection 
of our group-sized haloes (selected at redshift 1.0) with 'dark' 
substructure (crosses) and bright substructure (circles). The 5 
CLASS lenses found to have luminous satellite galaxies are plot- 
ted with solid circles; for B1127-I-385, the horizontal bar shows 
the range of separations when the lens redshift is varied from 0.5 
to 1. The histograms show the distribution of bright substructure 
found within the central 10 kpc (projected) as a function of 
separation and magnitude dilference. 



As haloes fall into a larger system, they are exposed to tidal 
forces and are stripped as they orbit the host system. The 
extent to which a halo is stripped d epends on resolutio n and 
the inner density profile of the halo (iMoore et al.lll996l ). The 
simulated subhaloes have artificially low density cores (due 
to force softening) that makes them more susceptible to tidal 
stripping. Including baryons (and gas cooling) will increase 
the central den sity and make the galaxy more resistan t to 
tidal stripping l|Moore et al.l ll996l : lMaccio et allbOOel ). al- 
though the coohng of baryons toward the central host will 
also increase the tidal forces experienced by subhaloes that 
come close to the centre. 

Nearly all of the satellite galaxies we find in the pro- 
jected central regions are tidally stripped 'orphan' galaxies 
(see Table |3}. The semi-analytic model we have used follows 
the orbit of galaxies which have lost their dark matter sub- 
halo, by assuming that they follow the motion of the most 
bound particle of the parent subhalo before it was destroyed 
(t his is shown to be a go od estimate of the subhalo's position 
bv lSpringel et af]|200ll ). Since the effects of dynamical fric- 
tion on the orbit of these galaxies is not considered in detail, 
ca ution is required w hen interpreting our results. As noted 
m ISales et al.1 (|2007l ). this may affect the overall abundance 
and radial distribution of these stripped haloes. 

Also, it is assumed that the 'orphan' galaxy remains 
completely undisturbed fo r a merging time, based on t he dy- 
namical friction formula of iBinnev fc Tremaind |l983) , until 
it merges with the central galaxy. This assumption may re- 
sult in an over-estimate of the number of 'orphan' ga laxies 
and their associated luminositv. lHenriques et al.l (|2003) take 
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the opposite approach and assume that all 'orphan' galaxies 
which have not merged with the central galaxy by 2: = 
are completely disrupted, and are responsible for the diffuse 
intracluster light. While their results suggest an improved 
match to the luminosity function in groups and clusters, the 
model is simplistic. They note that it is more likely that 
the disruption would happen gradually, and that the dense 
cores may survive for longer. The extent to which disruption 
would affect these 'orphan' galaxies remains unclear. How- 
ever, survival of these 'orp han' galaxie s (at l east to some 
extent) has been shown bv lWang et al.l (|2006l ) to be essen- 
tial in order to explain the observed correlation signal at 
small scales. This study used the Millennium Simulation to 
construct a new model of galaxy clustering. They found that 
if 'orphan' galaxies were excluded from the analysis, the cor- 
relation signal decreases at small scales in contrast to obser- 
vations. 

A related question is: if we were to increase the nu- 
merical resolution of the simulation, would the fraction of 
luminous satellites rise significantly? Clearly, the number 
of subhaloes (dark or luminous) must rise further since 
the subhalo mass function roughly fo llows a power law 
with dn/dM PC M~°, a = 1.7 - 1.9 (Moore ct al."l99d; 
Ghigna et al.ll2000l; [Pe Lucia et all [2004 : Gao ct al. 20049; 
Diemand et al.ll2007bl ). The lowest mass subhaloes we can 
resolve have circular velocities at the virial radius, Vc, of 
the order < 50 km s~^; haloes with Vc < 30 km s~^ may 
be inhibited from star formation b y the UV bac kground ra- 
diation (e.g. Rees, 1986 ; Efst athiou|[l992l : iThoul fc WeinberS 
ll996l : lGnedinl2000l V Thus many subhaloes may remain dark, 
and the fraction of bright subhaloes will not increase signifi- 
cantly. Increasing the resolution of the simulation would also 
mean that some of the 'orphan' galaxies would be resolved. 
Since only 'orphan' galaxies are allowed to merge with the 
central galaxy, increasing the resolution may prolong the 
lifetime of some of our 'orphan' galaxies. To fully under- 
stand the impact of this effect a more quantitative analysis 
is required. Clearly, a firm conclusion can only be reached 
with higher-resolution simulations with realistic treatment 
of the gas processes. 



4 SUMMARY AND DISCUSSION 

In summary, for the CLASS survey, approximately 5 of the 
22 primary lensing galaxies appear to have a faint com- 
panion within the projected central 5 kpc. The com- 
panions have luminosities of about 2 — 40% of the primary 
galaxy. We have studied host galaxies covering a compara- 
ble range of luminosities and host-to-satellite separations to 
the CLASS lenses, and found that the predicted fraction of 
galaxy- (group-) sized haloes hosting central luminous satel- 
lites (~ 3% (6%) at z = 0; ~ 11% (17%) at 2 = 1) is slightly 
lower than (but possibly consistent with) the observed value. 

While this fraction is largely independent of galaxy 
type, it is shown to increase with redshift. The Poisson prob- 
ability of detecting luminous substructure in 5 out of 22 
lenses, given that 3% of haloes host luminous substructure, 
is ~ 5 X 10~*. If 17% of haloes host luminous substructure, 
the probability of such a detection is ~ 0.14. Our prediction 
of the redshift and mass dependence appears to be roughly 
consistent with the data: three lenses with luminous satel- 



lites are in groups (see Section[3]2}, and all appear have red- 
shifts close to 1, higher than the median redshift (~ 0.6) of 
all CLASS lenses (see the right panel of Fig.|3J. One possibil- 
ity that we have not considered, is whether lensing galaxies 
are biased tracers of substructure; such bias may arise if sub- 
structure enhances the lensing cross-sections significantly. 
Previous studies, on cl uster scales, for gian t arcs indicates 
that the bias is small (|Hennawi et al.|[2007h : it remains to 
be seen whether this holds true for galaxy-scale lenses. Ob- 
servationally, the Sloan Lens ACS Survey (SLAGS) seems 
to indicate that the lensing galaxies at z ~ 0.2 are typical 
early- type galaxies (|Treu et al.ll2008l ). Another possibility is 
that some of the luminous 'satellites' are not associated with 
lensing galaxies at all, but just happen to be along the line 
of si ght dM ctcah 2005). 

I Shin fc Evans (2008, ') recently studied the effect of satel- 
lite galaxies on gravitational lensing flux ratios using an- 
alytic expressions for the host potential and the satellite 
galaxies. They use a spherically symmetric galaxy distribu- 
tion, and assume that the three-dimensional number density 
falls off like r~^'^ , comparable to the Milky Way. They show 
that the probability of a finding a large dwarf is about 10% 
within two Einstein radii and about 3% within one Einstein 
radius. We find that our z — results are consistent with 
this. 

The Millennium Simulation assumes a power-spectrum 
normalisation of erg = 0.9, slightly higher than the lat- 
est WMAP five-year result (|Komatsu et al.1 1200^ . where 
erg — 0.8. A lower value of erg will mean that haloes are 
expected to form later and be less concentrated. However, 
the impact of this parameter on our results is complicated. 
The semi-analytic models allow some fin e-tuning of parame - 
ters to match observations. For example. IWang et al.l (|2008t ) 
found no significant difference in the galaxy populations (at 
the redshift range relevant here) create d from semi-analytic 
models based on the W MAP one-year ( Spergel et al.l 120031 ) 
and WMAP three-year l|Spergel et al]|2007l ') erg values of 0.9 
and 0.722, provided suitable galaxy formation parameters 
were chosen (the difference becomes significant at high red- 
shift). We find that our results do not change significantly 
wh en based on the WMAP3 galaxy catalogue produced 
by I Wang et al.l (|2008l ) when the same merger timescale is 
adopted (as in their model C). However, in their model B 
(which has the same star fo rmation efficiency but a sh orter 
merger timescale than the iDe Lucia fc BlaizotI |2007| cata- 
logue) we find a factor of ~ 2 fewer haloes with central 
substructure. 

To summarise, while we find that the fraction of lu- 
minous satellites in group-sized haloes at z ~ 1 is roughly 
consistent with the observational data we caution that a firm 
conclusion can only be reached with higher-resolution sim- 
ulations involving a realistic treatment of the gas processes. 
At the same time, a larger sample of gravitational lenses will 
also be beneficial to constrain these models and allow more 
definitive conclusions on the properties of the substructure 
to be made. 
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